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We study the effect of four-particle ring exchange process on orbital antiferromagnetic state that 
occurs in some correlated electron systems in two dimensions. The primary question is whether the 
ring exchange process enhances or suppresses the orbital antiferromagnetic ordering. Using the fact 
that the orbital antiferromagnetic state arises in the large- TV limit of the SU(N) generalization of the 
t-J model, we consider the large- TV limit of the t-J-K model where K represents the four-particle 
, ring exchange term. The phase diagrams in the large-TV mean field theory are obtained for the 

half-filling and finite hole concentrations at zero temperature. It is found that the ring exchange in 
general favors dimerized states or the inhomogeneous orbital antiferromagnetic state, and suppresses 
the homogeneous orbital antiferromagnetic state. We compare our results with other related models 
of strongly correlated systems with ring exchange processes. 

PACS numbers: 71.10.Hf, 71.10.Fd, 71.30.+h, 75.10.Jm 
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Ring exchange interactions are often present in strongly correlated many-particle systems. In the Hubbard model 
of electrons with on-site repulsive interaction U , the large U limit leads to an antiferromagnetic exchange interaction 
J ~ t 2 /U between electrons at nearby sites where t is the hopping amplitude of electrons. It was shown that the 
next leading order term in the effective Hamiltonian is the four-particle ring exchange interaction with the strength 
£h ' K ~ t 3 /U 2 fjjj. The multi-particle ring exchange interactions have been studied extensively in solid 3 He 0. It was 
also suggested that the ring exchange interactions may play an important role in understanding the 'metal-insulator' 
transition in low density two-dimensional electron systems |3j . More recently, ring exchange models have been revisited 
in the context of the occurrence of novel quantum ground states in strongly correlated electrons, bosons, and spin 
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jhases where the elementary excitations carry fractions 



systems {§,||J||7) • In particular the possibility of fractionalized 

of quantum numbers of electrons/spins has been discussed |0,||. It was also suggested that in some soft-core boson 
models with ring exchange interactions, there may exist a novel phase which is neither a supcrfluid nor a conventional 
£SJ I insulator; the low energy excitations in this phase are described by a continuous set of wavevectors with zero energy, 
t-H ■ which may be called a Bose 'surface' in analogy to the familiar Fermi surface of electrons ||. On the other hand, 
a quantum Monte Carlo study of a quantum XY model with a ring exchange interaction (which may correspond to 
a hard core version of the boson model with a similar ring exchange considered in Ref. [g) did not find this novel 
phase, but found that new emerging ground state is a striped bond-plaquette ordered state at intermediate values of 
■ the ratio between the four-spin ring exchange and the nearest-neighbor Heisenberg exchange couplings |t]]. 

Independent of these developments, there have been great interest in finding novel broken symmetry phases in 
I ■ correlated electron systems. Notable examples are the proposals of various time-reversal symmetry broken phases for 
the explanation of the pseudogap phase of underdoped cuprates. HJ|,|lCj One of these examples involve short-ranged 
(fluctuating) Q or long-ranged Q orbital antiferromagnetic order; the long-range ordered state is often called the 
staggered flux phase [jll| or the d-density wave state S. This is a metallic state and characterized by alternating 
circulating currents in plaquettes of the square lattice while the charge density is uniform. Whether the short- 
. £^ | range or the long-range ordered orbital antiferromagnetism is responsible for the pseudogap in cuprates has been a 
subject of recent debates |l^ , |i3| ]. Mean field theory studies of some microscopic models have shown that the orbital 
antiferromagnetic state appears in some part of the phase diagram |^l| , ^4[ . Some ladder models also seem to allow 
the orbital antiferromagnetic phase (IB) . 

Ubiquitous presence of ring exchange interactions in strongly correlated systems leads to the question of the stability 
of the orbital antiferromagnetic state in the presence of ring exchange processes. This is a particularly interesting 
question given that the large U limit of the Hubbard model possesses the ring exchange term . Previous neutron 
scattering experiments on an undoped cuprate La 2 _ 2; Sr ;c Cu04 suggests that the presence of ring exchange is necessary 
to explain the spin wave spectra while the effect of ring exchange interaction may become weaker at finite 
concentration of mobile holes. In this paper, we investigate the effect of ring exchange process in a two-dimensional 
correlated electron system. In particular we would like to understand the effect of ring exchange process on the 
orbital antiferromagnetic state. A good starting point would be the large N limit of the SU(N) generalization of 
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the t-J model where the orbital antiferromagnetic state is already present as a possible ground state. The ground 
states discovered in the large N limit may be better understood as possible 'disordered' states of the zero temperature 
magnetically ordered phases in the physical limit N — 2 (such as spin antiferromagnetism) , as some parameter at 
T = changes or temperature is raised. The stability of different phases in the large N limit will be reflected in those 
situations as well. 

In this paper, we study the effect of the four-particle ring exchange interaction in the SU(N) generalization of the 
t-J-K model in the large N limit. Here K represents the strength of the four-particle ring exchange. The phase 
diagrams in the large N mean field theory are obtained at the half-filling and at finite hole concentration for varying 
strength of K and the Heisenberg exchange coupling J. We found that the ring exchange interaction in general 
suppresses the homogeneous orbital antiferromagnetic state and favors bond-centered dimerized phases with broken 
translational symmetry and/or the inhomogeneous orbital antiferromagnetic state |X7[ depending on the value of K. 

The rest of the paper is organized as follows. In section II, the SU(N) generalization of the t-J-K model and the 
corresponding large N mean field theory are discussed. The phase diagrams for various values of J and K at the 
half-filling and finite hole concentration are constructed in section III. Summary of the results, discussions on possible 
relations to other models with ring exchange, and conclusion can be found in section IV. 

II. THE MODEL 

We consider the t-J model with four-particle exchange couplings on an anisotropic triangular lattice with nearest- 
neighbor and next-nearest-neighbor interactions: 

H = -tx 2J c\ a °ia ~h C ]" ClQ ~ t3 X! C i Qc ka 

<ij> <<jl>> <<ik>> 

+ Ji ^2 (Si ■ Sj - ^rnn } ) + J 2 ^ (£j ■ Si - ^n 3 ni) + J 3 ^ 0i ■ S k - ^n ink ) 

<ij> ' <<jl>> ' <<ik>> 

+ ]T [^(^•^)(4-5,) + (5 i -5 1 )(5 j -S ? k )-A(S ? i -^k)(S ? j -S f 1 )], (1) 

<ijkl> 

where Ci a is the electron destruction operator of spin a =1, 1 on site i, Si = ^c\ a (7^Cif3 is the electron spin operator, 
and rii = c[ a Ci a (sum over repeated indices are assumed throughout this paper) is the electron number operator. 
Here J\ and t\ are the nearest-neighbor Heisenberg exchange coupling and hopping parameter while J2, J3 and £2,^3 
are the next-nearest-neighbor Heisenberg exchange coupling and hopping parameters respectively as shown in Figjl]. 
K and A are the strength of the four-particle exchange couplings, and the sum in the four-particle exchange is over 
all non-overlapping plaquettes labelled clockwise by the sites i, j, k and 1 (see Fig. 1). Note that at A = K. the 
four-particle exchange terms correspond to the conventional ring exchange process [3|. To get further insight, we 
consider more general four-particle exchange interactions by treating A as an independent free parameter [|18| . By 
doing so, we may at first investigate several simple limits where the model is easier to solve and the results in these 
limits can provide us useful guidance in studying the effect of various four-particle exchange couplings (including the 
ring exchange) on the orbital antiferromagnetic state. Then we will come back later to the case of more interesting 
and physical ring exchange coupling. 



FIG. 1. Anisotropic triangular lattice with three types of the bonds. 
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To make a progress, we apply a large- TV approach by generalizing the spin symmetry from SU(2) to SU(N) pl|19[. 
N = 2 corresponds to the physical SU(2) spin-1/2 system. The large- N approach has been applied to various strongly 
correlated systems |y],[ll],^(| . One advantage of this approach is that the mean field solutions become exact in the 
limit of N — > oo and thus it is possible to systematically calculate corrections in powers of 1/N. 

Let us briefly describe the SU(N) generalization of the Hamiltonian in Eq.|l|. In the SU(2) case, we can express 
the Hamiltonian in terms of fermionic operators by use of the following identity fLlfl : 

Si- Sj = ~^(4 a c ia )(cfc ip ) + ~rji - irnraj, (2) 

The SU (N) generalization of the SU (2) Hamiltonian is obtained by simply letting the spin index a run from 1 to N. 
We apply the slave-boson method to impose the constraint in the Hilbert space. In the slave-boson representation, the 
spins are represented by fermions f a , a = 1 • • • N which transform under the fundamental representation of SU(N); 
the holes are represented by spinless bosons b\ where c\ a = fi a b\. The local constraint is given by [^T| : 

fl a h a + b\h = ^. (3) 



In the N — > oo limit and at zero temperature, the average hole concentration 6 can be obtained from j2l| 

J_ .+„ . N 



-^E<// Q ^>=^( 1 -< 5 )> ( 4 ) 



where N s is the total number of sites in the system. In this paper, we only consider the case with uniform charge 
distribution which excludes the possibility of site-centered charge density wave or stripe state that was found in some 
related models . Investigation of the possibility of charge stripe phases is an interesting and yet a separate topic 
which is beyond the scope of this paper. Based on this assumption, the hole concentration is given by 5, and the 
electron density operator n- la is a constant c number, n\ = ^(1 — 5). Note that there is exactly one electron per site 
in the limit of the physical SU(2) spin-1/2 case at the half-filling. 

The SU(N) generalization of the Hamiltonian in Eq{| is therefore given by: 

H S u(N) = -hSj2fi a f^-t2S J2 ftf^-hS £ ft ha 

<ij> <<jl>> <<ik>> 

_ Z±_ I ft" f |2 _ J% \^ f \ot f \2_'Jz_ I f\ot f |2 

N l-'i •'J q| TV ^ j N 

<ij> <<jl>> <<ik>> 

+ E [^(l/i tQ! /j«l 2 [/^/^l 2 + l/^/icl 2 !//^/^! 2 ) - ^l/^/kcl 2 !/,^/^! 2 ], (5) 

<ijkl> 

where the Heisenberg exchange couplings Jy have been renormalized by the ring exchange interactions leading to the 
effective couplings: Jf = Ji(l - ^-(1 - <5) 2 ), Jf = J 2 (l + ^(1 - S) 2 ), and Jf = J 3 (l + ^(1 - S) 2 ). Note that we 
rescale both Jy and K by a factor of so that the Hamiltonian is of the order N . PH . Also, we have dropped the 
constant terms involving ni and nirij in the Hamiltonian. 

The above Hamiltonian contains eight- and four-fermion interactions. We may decouple these interactions into 
quadratic interactions by applying the Hubbard-Stratonovich transformations. We then perform the functional in- 
tegration over the Hubbard-Stratonovich fields in a saddle-point approximation The key point in the large- TV 
approach is that the saddle-point solution becomes exact in the N — > oo limit. The factorization of the spin-spin 
interactions can be done in two stages . We first break up the eight-fermion terms by introducing real Hubbard- 
Stratonovich fields $ij : 

^\th a ?\frh,\ 2 - iv[i$ ik % - ^i^/^i 2 - ^\fFhe\ 2 } (6) 

At the saddle-point $ij fields take the following values: 



- - — < |ft"f. |2 > 

$ik=^<|/ ; tQ /kar>, (7) 
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where $y and <J>ik represent the fields along the nearest-neighbor and next-nearest-neighbor sites respectively. 

After the first Hubbard-Stratonovich transformation, the Hamiltonian now still contains four-fermion interactions. 
We can decouple these terms by introducing the complex Hubbard-Stratonovich fields xy : 

-Jzii + ^)\f\ a h a \ 2 - (i + ^)(^# - § f} a fi» + 

-Jl(l + ^)\fl a h a \ 2 - (1 + ^)(^ " ^//"/k. + H.c). (8) 
where xy fields have the following saddle-point values : 

XH = § < /i ta /ja > • (9) 

After two Hubbard-Stratonovich transformations, the Hamiltonian is now quadratic in fermion operators. Note 
that at the saddle-point $y and xy fields have the following relations: 

2A 

^k^TT^IXikl 2 , (10) 

where ij and ik again are the nearest-neighbor and next-nearest-neighbor sites respectively. Therefore, we can express 
the mean-filed Hamiltonian purely in terms of x-fields. After integrating out the fermions, the resulting free energy 
is given by jll| : 
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E"MF r 6xY |2 | ,2 i i9| 1 2 \ 6-K . |2 , , 

~~/\T = 2^ [-TjSud^jl l^ kI l +M IXJkl ) + (je) 2 (je)2 IXik| |Xji| j 

<ijkl> ^ 1 ' \ 2/ \ 3/ 

\ - IXijl 2 \- Ixjil 2 \- IXikl 2 

je 2-~i je 2—1 je 

<ij> 1 <<jl>> 2 <<ik>> 3 

-^k{l + e[-M}, (11) 
' k 

where fj, is the chemical potential which fixes the number of the fermions, j3 = is the inverse of the temperature, 
Wk are the eigenvalues of the following Hamiltonian H\ : 

§- = -ti5j2fl a fi«-t2S £ fl a h a -h5 £ jf/ kQ 

<ij> <<jl>> <<ik>> 

+ E [(i-^?^)(-|// Q / ja +H.c.) 
<ijki> [J ^> iV 

i /i 4if|xil| Xjk ,+ a p tt \ 

+ i (Jf)3 ^~ aT^ /ka + H - c - } 
^ 2%ikp Xji tQ ,0-1 

+ (1 + ^Jf|J)(-^/i tQ /^ + *•«■)]• (12) 
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FIG. 2. The 2 x 2-site unit cell of the SU(N) t-J model on an anisotropic triangular lattice. Arrows indicate the orientation 
of the complex-valued Xij fields. 

Note that at the half-filling, the mean-field Hamiltonian is invariant under local U(l) gauge transformations as 
long as the x-fields transform as gauge fields @: xyO) -> e^'W-^^XijO). Away from half-filling, this local U(l) 
gauge symmetry is broken down to only the global U(l) gauge symmetry reflecting the conservation of total charge. 
The ground state of the system is found by solving the saddle-point equations, d ^ MF = 0, for Yij-fields. Extrema 

of Emf are found numerically with the simplex- annealing method p2|. To make progress in solving the model, we 
make the assumption that spontaneous symmetry breaking, if it occurs, does not lead to a unit cell larger than 4 
sites. Our choice of unit cell is shown in Fig.[| The 2x2 unit cell requires 16 different complex-valued x-fields (8 
Xi-fields on the square and 8 X2-fields along the diagonals). We work with a lattice of 40 x 40 sites and check that 
this is sufficiently large to accurately represent the thermodynamic limit. We check all the saddle-point solutions and 
make sure that the ground state is the one with the lowest free energy. 



III. ZERO TEMPERATURE PHASE DIAGRAM 



Possible ground states of the system may be classified into several different phases which are characterized by 
the symmetry of the xy- fields at the saddle-point. In the half-filled (5 = 0) case, it is important to classify the 
phases in a gauge-invariant fashion because there are many gauge-equivalent saddle-points. There are two important 
gauge-invariant quantities [^Tj| : 

(i) The squared amplitude |xij| 2 which is proportional to the spin-spin correlation function (Si • Sj). Modulations in 
|x| signal the presence of a bond-centered dimerization. 

(ii) The plaquette operator II = Xi2X23X34X4ij where 1, 2, 3, and 4 are sites on the corners of a unit plaquette. By 
identifying the phase of x as a spatial gauge field it is clear that the plaquette operator is gauge- invariant, and its 
phase measures the amount of magnetic flux penetrating the plaquette 11|. Different saddle points are therefore 
gauge equivalent if the plaquette operator has the same expectation value, even though the x-ficlds may be different. 
At half-filling the flux always equals to or ir (mod 27r), so one can always choose a gauge such that all the x-fields 
are purely real (ll) . 

Away from the 8 = limit, the saddle-point solutions are further classified in terms of whether or not they break 
time-reversal symmetry (T) where the plaquette operator has a nontrivial phase other than or 7r. The T-breaking 
state will generate real orbital currents going through the plaquettes in an alternating fashion [ fu|]2"c|l . Finally, as there 
are eight independent parameters (ti, £2, ^3 , Ji, ^2, ^3, K, and A) the phase diagram lives in a seven-dimensional 
space of their dimensionless ratios. Since we are interested in how the four-particle exchange couplings may affect 
the orbital antiferromagnetic state, we reduce the parameter space to a more manageable two-dimensional section 
by fixing the ratios of Jf/Jf) ^xl ^ti fa/fat W^j an d K/X such that the orbital antiferromagnetic state appears as 
one of the ground states of the system. Then by varying Jf/ti (or J£/t 2 ) and K/Jf (or A/Jf) we explore the phase 
diagram in a two-dimensional space. We summarize the phases which appear in the phase diagram in the following 
subsections. 
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A. Half-filling (S = 0) 



We find that in this limit the ground states are always dimerized in one way or another with spatial modulations 
in spin-spin correlations. Although the staggered flux phase (that becomes the orbital antiferromagnetic state at 
non-zero hole concentration) is a saddle-point solution, its free energy is always higher than the dimerized ground 
states. To sec this more clearly, we plot the free energy as a function of the four-particle exchange parameter K in the 
square lattice limit ( Jf = Jf = 0) with A = (see Fig](]). We compare the free energies of the saddle-point solutions 
which are competing with each other to be the true ground state. In this limit, we find the following competing 
phases: 

Dimer phase 

This is a fully dimerized phase which exhibits spin-Peierls order: = ^ 0, and all other ^-fields are zero; See 
Fig.|| for a sketch ^0|. It breaks the translational symmetry and have spatial modulations in spin-spin correlations. 
In fact, the system breaks up into decoupled dimerized spin chains. It is an insulating phase as there is a large gap 
in the energy spectrum at the Fermi energy. 



FIG. 3. The dimer phase with decoupled dimerized spin chains. Dark solid links indicate nonzero spin-spin correlations and 
dash links have zero spin-spin correlations. 



Box phase 

This insulating phase breaks the translational symmetry and consists of isolated plaquettes with enhanced spin-spin 
correlations pc| , p3| , p4| . The xi-helds are either or complex with \x\\ = \xl\ — IXil ~ Ixl I ^ 0, and \\\ \ — |Xil = 
Xil — IXil — 0. All xi fields are zero. See Fig.[|for a sketch. This phase is similar to the fully dimerized phase as 
the magnitudes of the xi-fields modulate on the lattice. The box phase does not break the time-reversal symmetry 
as the phase of the plaquette product X1X1X1X1 is either or it. There are no real orbital currents circulating since 
we can always choose a gauge to make x-helds real. 




FIG. 4. Box phase with or n flux around the plaquette. Dark solid links indicate nonzero spin-spin correlations and dash 
links have zero spin-spin correlations. The arrows indicate the complex-valued ^-fields. 



Staggered Flux Phase 

All %1-fields are equal, with an imaginary component in general, and all X2-fields along the diagonals are zero; see 
Fig.|| for a sketch. Like the box phase, the staggered flux phase at the half-filling does not break T-symmetry as the 
phase of the plaquette operator is either or tt. We can always make a gauge transformation such that all xi-ficlds 
are real. There are therefore no real orbital currents circulating around the plaquettes. This phase is semi- metallic as 
the density of states is small, and in fact vanishes linearly at the Fermi energy. The gauge fluctuations at sufficiently 
small- -/V are expected to drive the staggered flux phase into a Neel-ordered state p5| . 
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FIG. 5. Staggered flux phase in the square lattice limit with A = 0. The dark links on the square indicate strong spin-spin 
correlations and the dash links along the diagonals have zero spin-spin correlation. The arrows indicate the orientation of the 
complex- valued ^-fields. At half-filling, there is no real circulating orbital currents since the flux around the plaquette is either 
or 7T. Away from half-filling, however, the flux around the plaquette is neither nor it. There are therefore real orbital currents 
in this case circulating in an antiferromagnetic fashion as indicated by the arrows; thus it is an orbital antiferromagnet. 

We plot the free energy as a function of K/Jf as shown in Fig.|| We can see that the box and the dimer phases are 
the ground states for K < and K > respectively. The staggered flux phase is not the ground state in both cases. 
We have also checked the free energy as a function of A/J| in the other extreme limit where J-f = and K — (two 
inter-penetrating square lattice limit). We also find that the staggered flux phase is not the ground state for both 
A > and A < 0. 

-0.05 I 1 1 1 1 1 1 ' 1 




-0.2 1 1 1 1 1 1 1 1 1 

-2-10 1 2 

K / J" 

FIG. 6. Free energy of the half-filled t — J model on the square lattice as a function of Kj Jf for three competing saddle-point 
solutions. Note that the dimer and box phase are degenerate ground states at K — 0. 

From our results at the half-filling, it is clear that the four-particle exchange actually favors the dimerized phases 
over the staggered flux phase. The dimerized state which breaks translational symmetry is in fact a bond-centered 
charge density wave state fll}] . Recent Monte Carlo study of the quantum XY model with ring exchange indicates 
plaquette (stripe) order which is also a bond-centered charge density wave state as the ring exchange coupling exceeds 
the critical value 0). In fact, the box phase in our model in the square lattice limit with K < looks very similar 
to the plaquette phase in the quantum XY model mentioned above. Although there are differences between the XY 
model and our t-J model, these two models share the same structure of the spin-spin interactions and have similar 
ground states. Note that the four-particle exchange and the ring exchange couplings share the similar form of spin-spin 
interactions. We therefore expect that even though the staggered flux phase is not a ground state in our model at 
half-filling, our result at half-filling can still give us an important clue about the effects of the ring exchange couplings 
on the system at finite doping where the staggered flux phase does appear to be a ground state. Now we proceed to 
solve our model at finite doping. 
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B. Nonzero doping (5 = 0.05) 



Unlike the half-filled case where the staggered flux phase is never a ground state, at nonzero doping the hoping 
terms are expected to stabilize the staggered flux phase JTT) . We find the staggered flux phase is indeed the ground 
state of our system over a range of J\/t\ ratio. From our results at the half-filling, however, we know that the 
four-particle exchange favors dimerization over the staggered flux phase. We find it is also the case at finite doping. 
There is therefore the competition between the hopping terms and the four-particle exchange couplings which results 
in an interesting phase diagram. 

Again, we first work on the phase diagrams of two simple limits where the staggered flux phase appears as a stable 
ground state. We then move on to the more interesting ring exchange limit where all the couplings are nonzero. Here, 
we fix the doping concentration 5 to be 0.05. 



1. Square lattice limit with S = 0.05 

As shown in Fig.[l0| we find the following phases in our phase diagram: 
Dimer phase 

This dimer phase breaks translational symmetry and it exists at large Jf/t\ and positive K/Jf ratio. This phase 
is a partially dimerized phase as |x!| = |xfl > \xi\ = IXil > |x?| — Ixll > Ixll = IXil- AH X2-fields are zero. See Fig.0 
for a sketch. 



FIG. 7. The one-dimensional partially dimerized phase at S — 0.05. Dark solid links indicate strong spin-spin correlations, 
light solid and dark dash links indicate weak spin-spin correlations. The X2-fields on the diagonals are zero. 

Box phase 

This box phase next to the dimer phase also breaks the translational symmetry and exists at large Jf/ti ratio and at 
negative or small positive Kj Jf ratio. This phase is very similar to the box phase at half-filling except now it is partially 
dimerized. All xi-fields are complex and \x\\ = \xi\ — \Xi\ — Xi > Xi — \Xi\ — Ixll = IXil- The imaginary parts of 
Xi-fields have the similar relation: Im(x\) — Im(xi) = ^ TO (x?) — ^ m (Xi) > I m (Xi) — I m (Xi) — I m {Xi) = I m (Xi)- 
All X2-fields along the diagonals are zero. This box phase at finite doping breaks the time-reversal symmetry as the 
phase of the plaquette product is neither nor 7r in general. There are real orbital currents circulating around the 
plaquette. This corresponds to an inhomogeneous orbital antiferromagnetic state. See Fig.^j for a sketch. 




FIG. 8. Box phase on the square lattice at 5 = 0.05. The flux around the plaquette is neither nor ty. Dark solid links 
indicate larger magnitudes of the xi-helds or stronger spin-spin correlations and dark dash links indicate smaller magnitudes 
of the xi-helds or weaker spin-spin correlations. The arrows indicate the complex-valued xi-helds and also the directions of 
the orbital currents. The X2-fields along the diagonals are zero. We have checked that the orbital currents are conserved at 
each site. 



Staggered flux phase (Homogeneous orbital antiferromagnetic state) 

This phase exists at intermediate range of Jf/t\ ratio. It is very similar to the staggered flux phase at half- filling 
except that now it breaks T-symmetry as the phase of the plaquette operator is neither nor 7r; see Fig.O for a sketch. 
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There are therefore real orbital currents circulating around the plaquette. Note that unlike the box phase, the x-fields 
in the staggered flux phase have the same magnitude. 
Uniform phase 

This phase exists at small Jf/ti ratio where all xi-fields become negative real numbers and xi = Xi = " " = Xi> 
and all X2-fields are zero. See Fig.|| for a sketch. This phase preserves T-symmetry and since all x-fidds arc real, 
they simply renormalize the hopping parameter t\. The uniform phase is therefore a metallic Fermi liquid. Spin-spin 
correlations in the uniform phase decay with an inverse power law of the separation p0| . 



FIG. 9. Uniform phase in the square lattice limit. This phase does not have broken symmetries and it is a Fermi liquid. 



As shown in Fig|I(J, the general features of the phase diagram suggest that the region of the staggered flux phase 
(homogeneous orbital antiferromagnetic state) becomes narrower and that of the dimerized phase and/or the box 
phase (inhomogeneous orbital antiferromagnetic state) becomes wider as we increase both the positive and negative 
K/Jf ratios. This is consistent with what we find at the half-filled case where the four-particle exchange couplings 
enhance dimerization and suppress the staggered flux phase. Also, for a fixed Jf/ii ratio, the ground state may 
change from the box phase to the staggered flux phase and finally to the dimer phase as the ratio of K/ J\ is changed 
from negative to positive values. 




K / J° 

FIG. 10. Zero temperature phase diagram of the SU(N) t-J model in the square lattice limit as a function of Jf/ii and 
K/Jf. 



2. Two inter-penetrating square lattice limit with 8 — 0.05 



In this limit, the ground state phases are classified as follows: 
ID dimer phase 

This phase is a 'one-dimensional' partially dimerized phase which breaks up the system into nearly decoupled 
dimerized spin chains along one of the two diagonal directions: | X2 1 > Ixll > \xi\ = IX2 = Ixll > Ixll = Ixll = I X2 1 
(see Figjlin. This phase is the ground state for A < 0. Note that at half- filling, this phase is completely dimerized: 
1X2 1 = X2I 0j an d all other x-helds are zero. 
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FIG. 11. The 'one-dimensional' partially dimerized phase along the diagonals. Solid links indicate strong spin-spin correla- 
tions and dark dash links indicate weak spin-spin correlations. The xi-helds on the square are zero. 

2D dimer phase 

This phase is a 'two-dimensional' partially dimerized phase which breaks up the system into nearly decoupled 
dimerized spin chains along both of the diagonal directions: : X2I > \xl\ > Ixil = \xi\ = Ixil > X2 = X2I = IX2I 
(see Fig Jl^) . This phase is the ground state for A > 0. Note that at half- filling, this phase is completely dimerized: 
1X2 1 = X2I ^ 0' an d all other x-fields are zero. 




^ v 

FIG. 12. The 'two-dimensional' partially dimerized phase along the diagonals. Solid links indicate strong spin-spin correla- 
tions and dark dash links have weak spin-spin correlations. The xi-helds on the square are zero. 

Staggered flux phase (Homogeneous orbital antiferromagnetic phase) 

The staggered flux phase exists at intermediate values of J|/^2- All %2-fields are equal, with an imaginary component 
in general, and all xi-fields on the square are zero (see Fig.|l3|). This phase is very similar to the staggered flux phase 
on the square lattice except that now the flux goes through the plaquettes of two square lattices (rotated by 45°) 
inter-penetrating to each other. Like staggered flux phase on the square lattice at nonzero doping, this phase also 
breaks T-symmetry and there are real orbital currents running through the two inter-penetrating plaquettes as the 
phase of the plaquette operator is neither nor ir. 




FIG. 13. Staggered flux phase in the two inter-penetrating square lattice limit. The dark links along the diagonals indicate 
strong spin-spin correlations and the dash links on the square have zero spin-spin correlation. The arrows indicate the orientation 
of the complex- valued x-helds. The lattice can be also viewed in this limit as if there are two inter-penetrating square lattices 
rotated by 45° . At half-filling, there is no real circulating orbital currents since the flux around the plaquette of these rotated 
square lattices is either or n. Away from half-filling, however, the flux is neither nor n. Therefore, there are real orbital 
currents circulating in an antiferromagnetic fashion as indicated by the arrows. 

Uniform phase 

This is the ground state at small ratio. This metallic phase is similar to the uniform phase in the square 

lattice limit except that now \\ — ' ' ' = X2 7^ an d all xi-fields are zero (see Fig.[l4"|). 
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FIG. 14. Uniform phase in the two inter-penetrating square lattice limit at S — 0.05. This phase does not have broken 
symmetries and it is a Fermi liquid. 

We again plot the phase diagram as a function of Jf/i2 and A/Jf as shown in Fig.|l5|. The phase diagram suggests 
that the region of the staggered flux phase becomes narrower and that of the dimerized phases becomes wider as both 
the positive and negative A/ Jf ratios are increased. We reach the same conclusion as in the square lattice limit; the 
four-particle exchange couplings enhance dimerization and suppress the staggered flux phase (homogeneous orbital 
antifcrromagnetic state). Also, for a fixed Jf/ti ratio, there are quantum phase transitions between the ID dimer 
phase and the staggered flux phase as well as between staggered flux phase and the 2D dimer phase as the ratio of 
A/ Jf is varied. 

2 I . 1 . 1 . 1 . 1 




In addition to the above two extreme limits, we also look at the more interesting ring exchange limit where K = A 
and all Jy and iy are nonzero. We explore a wide range of the parameter space and find that the staggered flux 
phase only exists in the region close to the above two extreme limits. It is expected because the staggered flux phase 
is suppressed by geometrical frustration effect in other regions of the phase diagram |pof . We therefore restrict our 
attention to the parameter space close to these two limits where the staggered flux phase exists as a stable ground 
state. 

The phase diagrams of the system close to the square lattice limit and to the limit of two inter-penetrating square 
lattice are shown in Fig.[l6| and Fig]l7] respectively. They are very similar to the phase diagrams in their correspond- 
ing extreme limits except for the shifting of the phase boundaries. It is clear from our phase diagrams that the 
ring exchange coupling suppresses the homogeneous orbital antiferromagnetic state and favors the dimerized and/or 
inhomogeneous orbital antiferromagnetic state (box phase). 
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FIG. 16. Zero temperature phase diagram of the SU(N) t-J model at 8 — 0.05 in the ring exchange (K = A) limit. Here, 
J2/J1 = J3/J1 = 0.5 and t'2/ti — tz/t\ = 0.1. All the phases here are very similar to the corresponding phases in the square 
lattice limit except that now the X2-fields along the diagonals take small nonzero real values due to frustration. 




0.5 1 

K / i\ 

FIG. 17. Zero temperature phase diagram of the SU(N) t-J model at 8 — 0.05 in the ring exchange (K = A) limit. Here, 
J1/J2 = J1/J3 = 0.25 and ii/t2 = ti/t3 = 0.5. All the phases here are very similar to the corresponding phases in the limit 
of two inter-penetrating square lattice except that now the xi-ficlds on the square take small nonzero real values due to the 
nonzero couplings on the square links. 

IV. CONCLUSION 

In this paper, we studied the t-J model with four-particle exchange couplings given by Eq.l on an anisotropic 
triangular lattice as shown in Fig.l. Even though the usual ring exchange limit corresponds to K = A in Eq.l, we 
consider a more general model with K 7^ A in order to get further insights. The large- N mean field theory via the 
SU(N) generalization of the model is used to construct various phase diagrams. The mean field ground states in the 
A" — > 00 limit can be classified using the link variables Xij given by Eq.9. The relevant link variables are shown in 
Fig.2. 

We are mostly interested in the effect of the four-particle exchange couplings on the orbital antiferromagnetic state 
(staggered flux phase) that exists in some parts of the large parameter space of the model. It is found that it is useful 
to first consider two extreme limits; the square lattice limit (all the parameters describing the processes between the 



12 



next-nearest-neighbor sites are zero) and the limit of inter-penetrating square lattice (all the parameters describing 
the processes between the nearest-neighbor sites are zero). This is mainly because, even in the ring exchange limit 
K = A, the staggered flux phase exists as a stable ground states near these two limits; some small nonzero next- 
nearest-neighbor (nearest-neighbor) hopping and Heisenberg exchange couplings about the square lattice limit (the 
inter-penetrating square lattice limit). We studied these two extreme cases as well as the conventional ring exchange 
limit at half-filling and a finite hole concentration. 

Our results can be summarized as follows. In the half-filled case (S = 0), Fig. 6 shows the comparison of saddle-point 
energies of various phases in the square lattice limit. In this case, the staggered flux phase is never a ground state 
and cither the fully dimerized phase or the box phase is the ground state. One can see from Fig. 6 that the difference 
between the saddle-point energies of the staggered flux phase and the ground state (dimerized/box phase) becomes 
larger for large \K\ and/or |A|. At a finite hole concentration (5 = 0.05), the phase diagrams in the square lattice and 
the inter-penetrating square lattice limits are given by Fig. 10 and Fig. 15. Note that the finite hopping amplitudes 
can stabilize the staggered flux phase that has circulating orbital currents (orbital antiferromagnetic state) while the 
extremely large hopping amplitude eventually leads to the uniform phase. It is clear from these phase diagrams that 
the four-particle exchange processes favor the dimerized (these are 'partially' dimerized) phases or the inhomogeneous 
orbital antiferromagnetic state (box phase) over the homogeneous orbital antiferromagnetic state (staggered flux 
phase). This trend persists even in the conventional ring exchange limit K = A as can be seen from Fig. 16 and Fig. 17. 
Thus we conclude that the ring exchange or the four-particle exchange in general suppresses the staggered flux phase 
(homogeneous orbital antiferromagnetism) , and instead favors the dimerized phases and/or the inhomogeneous orbital 
antiferromagnetic state depending on the value of K. 

We note that a recent quantum Monte Carlo study j?J of a quantum XY model with a ring exchange process 
discovered a striped bond-plaquette ordered state at intermediate values of the ratio between the four-spin ring 
exchange and the nearest-neighbor Heisenberg exchange couplings. These results seem to be qualitatively consistent 
with our findings in the t-J-K model studied here even though there are obvious differences in two models. The 
existence of more exotic phases |||| like fractionalized phases in different parameter regimes of the t-J-K model 
would be an interesting subject of future study. 
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